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Introduction: 

Design  For  Blast  Mitigation 


MSTV 


MODELING  hnd  SIMULATION,  TESTING  rno  VALIDATION 


w 


Time  = 


Crew 


AlonBrill,  Boaz  Cohen  and  Paul  A.  Du  Bois,  SIMULATION 
AN  APC.  6th  European  LS-DYNA  Users'  Conference 


Composite  armor 
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Panel 


Landmine 


Restraint  system 
(seat  belt,  etc.) 


The  blast  mitigation  design  problem  can  be 

reduced  sub  problems  as  given 

Each  reduction  in  problem  formulation 

feeds  back  into  the  system  above 

Design  objectives  for  each  sub  problem  are 

selected  with  the  overall  problem  in  mind 

Vehicle 

-  Design  for  crew  and  critical 
component  survivability 

Sub  System 

-  Design  for  mechanical  isolation 
between  occupant  and  blast 

Component 

-  Design  for  minimum  energy  transfer 
from  blast  wave 

Sub  Component 

-  Design  for  energy  dissipation  and 
distribution 

Microstructure 

-  Define  damage  and  material 
parameters  for  energy  absorption 

OF  A  MINE  BLAST  EFFECT  ON  THE  OCCUPANTS  OF 
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Introduction: 
Injury  Criterion 


•  Injury  criteria  of  vehicle  occupants  due  to  mechanical  input  taken  as 
the  design  objective  of  the  vehicle  design  problem 

•  Blast  impulse  is  key  the  metric  which  drives  injury  occurrences 

•  Compressive  forces  and  vertical  acceleration  taken  to  be  defining  factor 
in  injury  accumulation 


HYBRID  III  Simulant 

Symbol  (units) 

Assessment  Reference  Values 

Response  Parameter 

Head  Injury  Criteria 

HIC 

750  ~5%  risk  of  brain  injury 

Head  resultant  acceleration 

A  (G) 

150  G  (2ms) 

Neck  forward  flexion  moment 

+  My  (N-m) 

190  N-m 

Neck  rearward  extension  moment 

-  My  (N-m) 

57  N-m 

Chest  resultant  acceleration 

A  (G) 

60  G  (3ms),  40  G  (7ms) 

Lumbar  spine  axial  compression  force 

Fz  (N) 

3800  N  (30ms),  6672  N  (0ms) 

Lumbar  spine  flexion  moment 

+  My  (N-m) 

1235  N-m 

Lumbar  spine  extension  moment 

-  My  (N-m) 

370  N-m 

Pelvis  vertical  acceleration 

Az  (G) 

15,  18,  23  G  (low,  med,  high  risk) 

Tibia  axial  compressive  force 

F(N) 

F/Fc  -  M/Mc  <  1 

combined  with  Tibia  bending  moment 

M  (N-m) 

where  Fc=35,584N  and  Mc=225N-m 

Femur  or  Tibia  axial  compression  force 

Fz  (N) 

7562  N  (10ms),  9074  N  (0ms) 

Occupant  Crash  Protection  Handbook  for  Tactical  Ground  Vehicles  2000 

Ala  Tabieiand  GauravNilakantan,  Reduction  of  Acceleration  Induced  Injuries  from  Mine  Blasts  under  Infantry  Vehicles 
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•  Topology  optimization  process  redistributes  material  in 
the  design  domain  to  obtain  a  concept  design 

•  Hybrid  Cellular  Automata  (HCA)  algorithm  using  uniform 
internal  energy  density  as  a  design  objective 

•  Nonlinear  transient  analysis,  utilizing  LS-Dyna  for  finite 
element  analysis  (FEA) 


mm 


min  ^  |  S,  -  S'  | 
s.t.  = 


<t<t. 


Topology  optimization  to  generate  concept  designs 
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•  A  continuum-based  topology  optimization 

-  First  utilized  for  bone  remodeling  (Tovar’04) 

-  Extend  bone  remodeling  technique  for  crashworthiness  design  ^107) 

•  HCA  =  Cellular  Automata  (CA)  +  FEM 

•  CAs  are  characterized  by  local  interactions  Neighborhoods 


Global  Formulation 

find  x 

s.t.  h(x)  =  0 

*(*)  <s 

H(x)  =  0 

C(x)  <  Q 

jc,-  e  { 0, 1 } ,  i  =  l,...,n, 


Local  Formulation 

find  x,- 

s.t.  >',(*,)  ->'*  =  0 

Xi  €  {0. 1 }, 


vonNeumann  (2D:  A/=4,  3D:  A/=6) 


Local  CA  rules  and  basic  control  theory 
is  used  to  distribute  material 


Moore  (2D:  N= 8,  3D:  A/=26) 
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HCA  Overview: 
Algorithm 


MODELING  hnd  SIMULATION,  TESTING  rno  VALIDATION 


-  Current  design 

> 

$*(0) 

f(*Q) 

Apply  HCA  material 
update  rules 

^  xfk\  S W 

^Convergenc^\ 

Jpdate  globa 
setpoint 

no  A 

yes 

-H  New  design 


Mass  control 


Initial  design 


p(0)=x(0|p0 


CT7(°)=*(0)ayo  E^=x^E, 


m=xfP)E, 

'  I 
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Crash  or  Blast 
analysis 


U(xP>) 


Update  material 
distribution 
using  HCA  rule 

Ax  =  /(U,U*) 


X(k+ 1) 


*(*-!) po  \ 

(  £(*+!)=  xf.**l)E0  \ 

\  Er"-ot"^hj 

New  material  model 
'Uniform  Internal  Energy  Density" 
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•  Field  Variable: 

-  Original  crasHCA  algorithm  only  utilized  Internal  Energy  (IE)  at  the 
final  time  step 

•  IE  at  the  final  time  is  highly  dependent  on  the  simulation  termination  time 

•  Resulting  topology  is  drastically  different  depending  on  the  selected  end  time 

-  Changed  method  for  blast  to  use  the  IE  at  all  time  steps. 


-  Will  utilize  the  concept  of  a  fully  stressed  design  as  implemented  in 
the  Crash  version  of  the  HCA  algorithm.  m' 
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jjjj  ^  Modification  of  HCA  for  Blast:  MSTV  f 

Johnson-Cook  Material  Model  MODELING  hno SIMULATION,  TESTING  rnd  VALIDATION 


Material  Card  Selection 

•  Piecewise-linear  elastic 
plastic  material  card: 

-  Quasi-static 

-  Hardening 

-  Plastic  deformation 

•  Johnson-Cook: 

-  Can  be  used  for  dynamic 
loading  situations 

-  Strain  rate  effects 

-  Temperature  effects 

E  =  E0xp  and  G  =  G0x" 

C7  =  [A  +  f?£"][l  +  Clne][l  -T  *"'] 
==  A=A0x‘>,  B  =  Sq-vT  C  =  Co x4 


Johnson-Cook:  Effect  of 
density  on  Internal  Energy 


0  - i - 1 - 1 - i - i 

0  0.2  0.4  0.6  0.8  1 

Relative  Density  (x) 

- GVSETS 
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Modification  of  HCAfor  Blast: 
CONWEP  Blast  Model 
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msp; 


MODELING  hnd  SIMULATION,  TESTING  rno  VALIDATION 


Load  Type: 

-  Began  using  the  CONWEP  algorithm  for  the  blast  model  in 
the  3-D  solid  element  HCA  method. 


•  Quick  Analysis  time  (relative  to  MMALE) 

•  Required  minimal  changes  to  the  HCA  algorithm 

-  The  objective  is  to  design  substructure  that  responds  to  a 
blast  event  in  a  desired  manner.  CONWEP  can  be  used  in 
this  scenario  since  we  are  only  looking  at  the  response  of  a 
small  piece  of  structure  rather  than  the  whole  object. 


S 


Az 
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Figure  1 .  Definition  of  variables  in  the  US  Army  TACOM  Impulse  Model 
(Adapted  from  Westine  et  a!.,  1985). 


UNCLASSIFIED 


GVSETS 


Implementation 


MSTV 


MODELING  mid  SIMULATION,  TESTING  rnd  VALIDATION 


m 


As  a  proof  of  concept,  a  rectangular  design  domain  was  created  to  represent  a 
piece  of  armor. 


-  Design  domain  is  40  x  40  x  10  cm 
aluminum  (represents  armor 
substructure) 

-  Top  layer  is  40  x  40  x  1  cm  ceramic 
(represents  ceramic  top  plate) 

-  Domain  and  top  plate  have  fixed  x,  y,  and 
z  displacement  boundary  conditions  on 
all  sides. 

-  Blast  is  positioned  100  cm  up  from  origin 
(89  cm  from  top  center  of  plate) 

-  Hourglass  control  is  included  to  help 
prevent  complex  sound  speeds  arising  in 
low  density  elements 

-  The  target  mass  is  set  to  be  50%  of  a  full 
design  domain 


Aluminum  Design  Domain 


•  Generated  Topology  to  be  compared  against 
a  baseline  model  that  is  full  density,  but  half 
as  thick. 

•  The  top  of  the  baseline  design  will  be  94  cm 
from  the  blast  source  (i.e.  the  base  of  the 
domain  will  be  the  same  distance  in  both 
cases) 


UNCLASSIFIED 


GVSETS 


'SKUffiWfift' 


Results: 

Integrated  IE  Objective 


MSTV 

MODELING  hnd  SIMULATION,  TESTING  rnd  VALIDATION 


Model  mesh  size  to  0.5  cm  (a  symmetric  quarter  of  the  domain  was  run) 
Simulation  end  time  800  |js 
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Results: 

MSTV  | 

r*i 

Integrated  IE  Objective 

MODELING  hnd  SIMULATION,  TESTING  rno  VALIDATION  '  ' 

0 

y  y 

Acceleration  Profile  of  iteration  4\j|\|CLASSIFIEfr:ce'eratlon  Pro^e  of  iteration  30 


Resulting  topology  has  mass  where  it  would  be  expected  and 
satisfies  the  mass  target  constraint. 

There  is  an  order  of  magnitude  improvement  in  the  mean  nodal 
acceleration  of  the  bottom  of  the  design  domain  versus  the  baseline 
case. 

Peak  acceleration  is  misleading  because  of  nodes  that  are  being 
unrealistically  accelerated  relative  to  their  neighbors 


HCA  Topology 
Baseline 


20  30 

iteration 


HCA  Topology 
Baseline 


0  25 


IJ  I  5 


X:43 


0.05 


20  30 

iteration 
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•  This  investigation  showed  that  the  HCA  algorithm  could 
be  modified  to  produce  topologies  that  help  to  mitigate 
the  acceleration  transferred  to  the  occupant  from  a  blast 
loading 

•  Future  work: 

-  Investigate  further  the  use  of  IE  as  the  field  variable  in  the 
optimization  process 

-  Investigation  of  other  field  variables  to  drive  the  optimization  that 
are  more  appropriately  related  to  acceleration 

-  Mesh  refinement  study 

-  Continued  work  to  improve  convergence  and  to  mitigate  errors  in 
the  LS-DYNA  runs  (i.e.  complex  sound  speeds  arising  in  low 
density  elements) 


TOrararrK  et.  an  nrjrorew 

Cellular  AutUrlatQi-A^I^j7]l/fe0c/7an/ca/  Design,  2009  (in  review). 
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Verification  of  Monotonic  Relationship 
between  SED  and  Mass  Density 

i'if,  j 


MODELING  hnd  SIMULATION,  TESTING  rnd  VALIDATION 


i  ne  necewise-unear  eiasnc=piasnc  moaei  was  snown  DyxrrrKaiei 
to  have  a  monotonic  relationship  between  SED  and  mass  density 
under  the  SIMP  penalization  method 

A  similar  study  was  conducted  to  determine  if  penalizing  the 
Johnson-Cook  model  also  yielded  a  monotonic  relationship  between 
SED  and  mass  density. 

-  Setup  as  a  single  solid  LS-DYNA  cube  element  under  a  rapid 
fixed  loading 

As  in  the  standard  SIMP  scheme,  elastic  modulus  [and  shear  modulus)  is 

-  Mass  ancrpenalization  factors  are  varied. 
penalizecTaccoraing  Fo: 


E  =  E0xp  and  G  =  G0xp 


Johnson-Cook  model  calculates  a  von-mises  flow  stress  according  to: 

a  =  [A  +  B£n][l  +  Clne][l-T  *m ] 

Penalizing  this  von-mises  flow  stress  is  akin  to  penalizing  the  yield  stress.  This 
is  done  by  penalizing  the  parameters  A,  B, and  C. 


3atel,  N.MV  B.S.  Kang,  J.E.  Renaud&  A.  Tovar.  Crashworthiness  design  using  topology  °f  Mechanical  Design.  Accepted  for  publication:  December  2008 
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Johnson-Cook  Material  Model:  Penalization 
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•  As  in  the  standard  SIMP  scheme,  elastic  modulus  (and 
shear  modulus)  is  penalized  according  to: 

E  =  E0xp  and  G  =  G0xp 

•  The  Johnson-Cook  model  calculates  a  von-mises  flow 
stress  according  to. 

a  =  [  A  +  Be"  ] [1  +  C  In e] [1  -  T  ] 

•  Penalizing  this  von-mises  flow  stress  is  akin  to 
penalizing  the  yield  stress.  This  is  done  by  penalizing 
the  parameters  A,  B, and  C. 
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Internal  Energy  Internal  Energy 


ohnson-Cook  Material  Model:  Fixed  Load 


•  Under  a  constant  load,  the  material  does  not  behave  monotonically  under  any 
penalization  scheme 

•  Depending  on  choice  of  p  and  q,  we  may  have  to  significantly  increase  the  minimum 
WjSntjn^ty  aH°wed  in  the  CA  and  FE  models 
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Johnson-Cook  Material  Model:  Fixed  Displacement 


i 


A 


q  = 


x  10 


Relative  Density  (x) 


q=  0.7 


x  10 


Relative  Density  (x) 


A 


q=  0.6 


x  10 


Relative  Density  (x) 


A 


•  Under  fixed  displacement  the  IED  appears  to  have  a  monotonic  relationship  with  relative 
density 

|||  Blast  loadingsjnowever,  are  not  fixed  displacement  problems. _ 

HBHIz  UNCLASSIFIED 
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Modification  of  HCAfor  Blast 
CONWEP  Blast  Model 


2' 

f 

F/ 


•  CONWEP  blast  model:  Load-Blast  function 
in  Ls-Dyna,  is  an  implementation  of  the 
hemispherical  blast  models  of  Kingery  and 
Bulmash. 

•  Empirical  blast-loading  model  rather  than 
explicitly  simulating  the  progress  of  the 
shock  wave  from  the  high  explosive 
through  the  air  and  its  interaction  with  the 
structure 

•  Does  not  account  for  pressure  confinement 
properties  provided  by  imbedding 
explosive  charge  in  soil. 

•  Scaling  charge  sizes  for  better  agreement 
accepted,  but  applications  to  complete 
structures  limited  due  to  improperly 
modeled  load  distributions 

•  More  complex  structures  and  interaction  of 
detonation  products  and  debris  requires  a 
more  sophisticated  fluid  structure 
formulation. 

Victnjjn  Chgite-- _ 

Kingery,  C.  N.,  and  G.  Bulmash.  "Air-Blast  Parameters  from  TNT  Spherical  Air  Burst  and 

Hemispherical  Surface  ui  TJmTy*!  a  1 1  istiHtesearcT^!lDorato^ 

A^^mWW^Ground,  MD,  April  1984 
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MODELING  hnd  SIMULATION,  TESTING  rnd  VALIDATION 


P{r)  =  Pt  •  cos2  6  +  Pt  •  (1  +  cos2  9-2 cos 6) 


menRnn^fRTTy7Tun3amentaI^fffotectiv^esigfT 

55-1,  Washington,  D.C.,  November  1986. 
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•  Introduction 

•  Overview  of  Hybrid  Cellular  Automata  (HCA) 

•  Methodology 

-  Field  Variable 

-  Material  Model 

-  Blast  Model 

•  Implementation 

•  Results 
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